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Abstract 

In this paper we study a non-equilibrium resonant tunnelling problem where a non-interacting 
quantum dot is connected to two leads, one being the edge of an interacting 2-D topological 
insulator (Luttinger liquid) and the other being a usual Fermi liquid. We show that the current 
passing through the system can be expressed in terms of a non-equilibrium local single-particle 
Green's function of the Luttinger liquid lead which can be analysed using standard bosonization- 
Renormalization Group (RG) technique. In particular, some exact results can be extracted in 
the small bias limit with repulsive electron-electron interaction. A simple formula which captures 
the qualitative feature of the I-V relation over whole temperature and voltage bias range is being 
proposed and studied. 

PACS numbers: 72.10.Bg, 73.40.Gk, 73.63.Kv, 71.10.Pm 
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Figure 1: The configuration 

The discovery of topological insulators [U [2] has generated much interests and activities in 
the condensed matter physics community. In the case of 2-D topological insulators (Quan- 
tum Spin-Hall systems), helical edge states respecting time-reversal symmetry [3] exist and 
provide an example of interacting Luttinger liquid which can be studied experimentally up 
to room temperature. The transport properties of interacting one-dimensional systems have 
been studied extensively in the literature[U-[6] , especially after the discovery of edge states 
in fractional quantum hall (FQH) liquids [7j (chiral Luttinger liquids). For example, the 
problem of non-interacting quantum dot connected to leads of FQH edge states was studied 
by Chamon and WenjB], and later by Furusaki [9] at temperatures higher than the tunnelling 
strength. 

The emergence of topological insulator makes these 1-D theoretical models realizable. In 
this paper we study a simple model of a non-interacting quantum dot connected to two leads, 
one being the edge of 2-D topological insulator with electron-electron interaction, the other 
being a normal Fermi liquid lead (Fig{T]). We show here that the non-linear I-V relation 
of the system can be expressed in terms of a non-equilibrium local single-particle Green's 
function of the topological insulator lead which can be analyzed in the limit of small voltage 
bias and low temperature using a conventional renormalization-group (RG) approach. A 
few exact results are extracted from the RG analysis and a simple formula which captures 
the qualitative feature of the I-V relation over whole temperature and voltage bias range is 
being proposed and studied in this paper. 



2 



We start with the model Hamiltonian H = Hqi + Hq 2 + Hod + Ht, where 



H i = -iv F J dx(ip{ Rt1 .(x)d x ipi R ^(x) 
+ f/^W*>W*) 

+^(x)^(x)) 2 (la) 

and 

H02 = 2J £2fecJ fco .c 2fc)Cr (lb) 

k.a 

are terms in the Hamiltonian describing the edge of a 2-D topological insulator (H i) and 
the Fermi liquid lead (H 02 ), respectively. 

Hod = eo 5^ <%d a (lc) 
(j 

is the Hamiltonian describing the quantum dot and 

H T = ^{^(0)^ + ^(0)^ + ^.0} 

+T 2 {4+.(0)rf t + + M (Id) 

describes electron tunnelling between the dot and the leads. We note that since the prop- 
agation and spin directions of electrons in the topological insulator lead are tied together 
(helical edge), we can forget its spin suffix and the helical edges states of topological insulator 
behave like a spinless Luttinger Liquid [TUl [11] . 

H i can be transformed to a free boson system by standard bosonization technique 



where the fermion correlation function g\(t, t') = —i ( Tc^i iO -(0, t)ipf^(0, t') ) can be com- 

\ / Hoi 



puted straightforwardly |14j . In bosonization theory a Luttinger liquid is usually charac- 
terized by the interaction strength parameter K = (1 + ^;) _1 ^ 2 ; here we further define 
fl = \{K + — 2) which is a parameter we shall use frequently in the following. 

The DC current / flowing from lead 1 to 2 in the above system can be expressed as[ 

'dN 2a 



y v 

^ \ dt 

-ieT 2 J2(ti,Mt)l(t) -d+(t)Vv(0,t)) (2) 
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where I2 is the current flowing from lead 2 to the quantum dot. Since lead 2 is non- 
interacting, we can eliminate the lead-2 electron operators through their equation of motion. 
Following Ref.[IS] we obtain after some straightforward algebra 

/ = -e\T 2 \ 2 J2 I dep 2 (e)[2f(e-fi 2 )Im{G^ a (e)} 

a J 

+Im{G^ r ( e )}] (3) 

where G dda (e) and G dda (e) are the Fourier Transform of the "retarded" and "less" com- 
ponents of the on-site (Keldysh) Green's function of the quantum dot, Gdd,a(t,t') = 
—i(Tcd a (f)d^(t'))i where t and t' are time-ordered along a closed time contour from 
t = —00 to t = 00 and then back to t = —00. /^(s) and \x 2 are the density of states 
(which will be taken to be a constant later) and chemical potential of lead 2, respectively. 
f(e) is the Fermi-distribution function. Notice that the current is expressed solely in terms 
of a (non-equilibrium) local single-particle Green's function of the system. This result is 
possible because lead 2 is non-interacting. 

To calculate Gdd,a we follow the Keldysh path integral formalism [12] . The action of the 
system is [16] 

S = J dt[J dx y^{^i ;CT (x, t)i(d t + avFd x )ip 1} a(x, t)} 
dx(^2 VvOe, t)ipi )t r(x, t)) 2 

a 

+ ^{Tn/v(0, t)d a (t) + T 2 ^ 2 , CT (0, t)d a (t) + c.c} 

(7 

+ ^d a (t)(id t -£ )d a (t)+J2{Ut)d*(t) + c.c}\ 

a a 

(4) 

where a = ±1 and we have suppressed the propagation direction indices of the helical 
edge states for brevity. {£ a , are external source fields introduced to generate Gdd,a{t,t'). 
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Integrating out the fermionic fields {c 2 k, a ,c 2 k,a} and {d a ,d a } we obtain 
S e ff = J dt[J dx y^{^i ;CT (x, t)i(d t + avFd x )ipi t(T (x, t)} 

-| J dxC^'lp!^,^!^,^) 2 ] 

[^(t 2 )+T 1 ^ lj(T (0,t 2 )] (5) 

which is an effective action for lead-1 electrons only. gDa{t, t') is the Green's function of the 
quantum dot evaluated under the "dot+lead 2" action with 

9nl A (") = ,. * (6) 



w - e ± «r 2 

2iT 2 f(u - /i 2 ) 



( w _ eo )2 + r 2. 

where ri( 2 ) = 7rpi(2)|2~i(2)| 2 is the tunnelling width from the dot to lead 1(2). G^ a {t y t') can 
be obtained by taking functional derivatives of the generating functional with respect to the 
{iaiio} fields; we obtain 

G Mj a(t,t') = g D a(t,t') + |Ti| 2 / dt 1 dt2gD (r (t,ti) x 

Jc 

9laeff(t 1 ,t 2 )gDa(t 2 ,t') (7) 

where 

gi aeff (t,t') = -z (^^(0,^(0, f 



e// 

is the Green's function of lead-1 electrons evaluated at x = according to the effective 
action ^ in the absence of the external source terms. Different Keldysh components of 
Gdd,a can be extracted from Q using Langreth's sum rules [T7j. Combining Eqs. ^ and 
(J7j), we obtain for the tunnelling current ([3]), 



{Im^< e// H + 2/(w - ^ 2 )Ims* (u;)}. (8) 



^o) 2 + r| 

where the only unknown is the effective green's function gi a eff(t, t'). 

To evaluate gi ae ff(t,t') we express it in the form of a standard Dyson's equation 

9laeff{t,t') = gia{t,t')+ / dtldt 2 g la (t, tl)E ff ,T(tl, ^) 

(*2,0 (9) 



where gi a (t,t') is the Green's function of lead-1 electrons at x = evaluated at T\ = 0, and 
is the standard x = Green's function of spinless Luttinger liquid with interaction strength 
The self-energy term in ^ represents correction coming from the tunnelling term 
(last term in Eq.(5)). The tunnelling current ^ can be expressed in terms of the Fourier- 
transformed self-energy Eo-.-rfw) as 



l^iH^I 2 ^ J dcop 2 x 



2(/(w - AUi) - /(w - A*i))«£(w) - |<?f» f (£j r M + 2/(w - /ia)Sj T (w)) Y1 
Imj ^ '■ [10 

|i-<Ms^M| 2 ((c-e ) 2 + rl) 

We observe that the tunnelling current is completely determined by the self-energy 
function £ CTj y(u;). In the absence of electron-electron interaction £ CT) y(u;) = £^°y(u;) = 
|T 1 | 2 5f£ )cr (u;) and our main job here is to understand how E^(o;) is renormalized by the 
electron-electron interaction. 

In the following we shall study E CTj r in the small bias, low temperature limit T, /^l << 
Ti(2) using a standard bosonization-RG analysis jl]. In this limit we may keep only the long- 
time behavior of gDa{t) in Eq. ^ and forget about the more complicated intermediate time 
behaviors, i.e. we approximate 

g R(A) {t) = e -ot-r 2 |*| Uoo ^ 0; (n) 
T 2 e^ 2 * 



9dM ~ ^ 



To proceed further we compare the present problem with the problem of directly tun- 
nelling between a Luttinger liquid and a Fermi liquid through a simple tunnelling junction 
barrier t. The Fermi liquid fields can be integrated out as what we have done to derive 
Eq. The only difference in the direct tunnelling problem is that gDa{t) in Eq. ^ is 
replaced by g$ A \t) ~ =FwriV(0)5(t) and g^(t) ~ 2ziV(0)e^ 2t /i where JV(0) is the density of 



states on the Fermi surface. Comparing with Eq. (11) We see that the resonant tunnelling 



problem reduces to the direct tunnelling problem in this limit if we replace the dimensionless 



tunnelling parameter tN(0) — > £ ^ + r 2 anc ^ the sca hng behavior of the self-energy in the 
present problem can be inferred from the corresponding direct tunnelling problem which has 
been analyzed using well-developed bosonization-RG technique. We obtain immediately 



T X {E) ~ (E/E F ) 



l-K, 



eff 



in this approximation, where K e ff = 2K/(1 + -^)[9J- In particular for repulsive interaction 
(K < 1), Ti(E — y 0) —7-0 and the self-energy scales to zero in the infrared regime, which 
makes perturbative RG applicable. In this case it is sufficient to approximate S CT) r ~ 
with T\ — > T\(E). Physically, the vanishing of self-energy is an alternative way to express 
the well known result that the tunnelling between the Luttinger liquid lead and the Fermi 
liquid lead vanishes in the infrared limit [H [9]. We emphasize here that the renormalization 
of T\ is restricted only to the self-energy term and does not appear in other places in the 
current expression. With this approximation we obtain 



-2e 



7T 



where E = max{T, \oj— Eq. (12) is expected to be reliable in both the low temperature, 
small bias limit (T, — /x 2 |) << ^1(2) and in the high temperature, large bias limit (T, \fii — 
/X2I) ~ Ep where T\(E) — > T\ and the renormalization of T% becomes unimportant. We shall 
analyze the current using this approximate formula in the following. 

To study the tunnelling characteristics, we define eV = Hi — [i\ and e = Eq — \i\ and look 
at the differential conductance as function of V at different temperature regimes. Using 
the result that gi a (oS) ~ E e and approximating the Lorentzian function by a 5-function, we 
find that at very high-temperature regime (T ^> ri )2 ,£o), T±(E) — > T\ and J£ scales with 
temperature as T l3 ~ 1 , in agreement with Furusaki's result j9]. On the other hand, the low 
temperature regime T r 12 can be sub-divided into two regions: when eV T, the linear 
differential conductance scales as T e (l + aT 7 ), where 7 = min{l, $ + 2(1 — K e ff)} and a is 
a temperature-independent constant; but when eV ^> T — > 0, the current-voltage relation 
becomes non- linear and becomes V— dependent (see discussion below) at small V. The 
different scaling behaviors are summarized in Figj2} 

Next we consider (non-linear) differential conductance at zero temperature. First, we 
observe that in the off-resonance region (|eV|,r 1)2 e), the differential conductance goes 
as 1^1^(1 + a'lV^I 7 ), where a' (different from a) is a temperature-independent constant. The 
first term ( | V" | 6> ) is in agreement with previous result [8] and the second term |^/| 6I+7 is the 
leading correction coming from scaling of T\. Our RG amalysis suggests that the leading 
order {\V\ e and T e ) scalings are exact as long as g^.(a;)S^ T (a;)| £i ,_ > o — > which is the case 
for repulsive electron-electron interaction and is independent of the detailed structure of 




Figure 2: The leading-order temperature dependence of differential conductance in different 
regimes. Regime 1: T >> Ti ^; Regime 2: \eV\ « T « Ti^; Regime 3: T « \eV\. 

dV 



Vie 1 

Figure 3: (color online) Differential conductance at zero temperature with e = 

E^ T (tu). The self-energy gives rise only to higher-order corrections to V/T scalings at low 
temperature and small voltage bias. 

We also observe a resonance peak of differential conductance located around eV ~ e and 
the differential conductance at V — is always zero as long as lead 1 is interacting. This is 
a character of Luttinger liquids which is very different for non-interacting electrons (K = 1), 
where ^\v^o approaches a constant. Finally, when V is far away from the resonant point, 
the differential conductance falls as |y|^-2xmax{i^+2(i-x e// )}_ 

To conclude, we have extended Meir-Wingreen's current formula to the case with one 
interacting lead in this paper, and have used it to calculate the tunnelling current from the 
edge of a 2D topological insulator through a quantum dot to a normal Fermi liquid lead. The 
formulation allows us to construct a current expression in terms of the self-energy of a local 
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Green's function. The differential conductance in different temperature regimes is analyzed 
using perturbative RG in this paper, which is believed to be reliable when both temperature 
and voltage bias are much smaller than IV2) and if the electron-electron interaction in lead 
f is repulsive. Based on the RG result, an approximate formula for the current qualitatively 
valid over whole temperature/volatge range is proposed. The formula reproduces results 
of earlier works at high temperature and produces exact results at low temperature for 
repulsive electron-electron interaction and small bias. For attractive interaction the scaling 
breaks down at low enough energy E suggesting that qualitatively new dl/dV behavior 
is expected at low energy. Our approach offers a new theoretical tool of analysing (non- 
equilibrium) DC transports which can be extended to other systems with both Luttinger 
and Fermi liquid leads. 

We thank K. T. Law, Zhengxin Liu and C. Chan for useful discussions. This work is 
supported by HKRGC grant HKUST3/CRF/09. 
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